3.2.48 \(\int (c+d x+e x^2+f x^3) (a+b x^4)^3 \, dx\) [148]

3.2.48.1 Optimal result
3.2.48.2 Mathematica [A] (verified)
3.2.48.3 Rubi [A] (verified)
3.2.48.4 Maple [A] (verified)
3.2.48.5 Fricas [A] (verification not implemented)
3.2.48.6 Sympy [A] (verification not implemented)
3.2.48.7 Maxima [A] (verification not implemented)
3.2.48.8 Giac [A] (verification not implemented)
3.2.48.9 Mupad [B] (verification not implemented)

3.2.48.1 Optimal result

Integrand size = 25, antiderivative size = 151 \[ \int \left (c+d x+e x^2+f x^3\right ) \left (a+b x^4\right )^3 \, dx=a^3 c x+\frac {1}{2} a^3 d x^2+\frac {1}{3} a^3 e x^3+\frac {3}{5} a^2 b c x^5+\frac {1}{2} a^2 b d x^6+\frac {3}{7} a^2 b e x^7+\frac {1}{3} a b^2 c x^9+\frac {3}{10} a b^2 d x^{10}+\frac {3}{11} a b^2 e x^{11}+\frac {1}{13} b^3 c x^{13}+\frac {1}{14} b^3 d x^{14}+\frac {1}{15} b^3 e x^{15}+\frac {f \left (a+b x^4\right )^4}{16 b} \]

output
a^3*c*x+1/2*a^3*d*x^2+1/3*a^3*e*x^3+3/5*a^2*b*c*x^5+1/2*a^2*b*d*x^6+3/7*a^ 
2*b*e*x^7+1/3*a*b^2*c*x^9+3/10*a*b^2*d*x^10+3/11*a*b^2*e*x^11+1/13*b^3*c*x 
^13+1/14*b^3*d*x^14+1/15*b^3*e*x^15+1/16*f*(b*x^4+a)^4/b
 
3.2.48.2 Mathematica [A] (verified)

Time = 0.01 (sec) , antiderivative size = 180, normalized size of antiderivative = 1.19 \[ \int \left (c+d x+e x^2+f x^3\right ) \left (a+b x^4\right )^3 \, dx=a^3 c x+\frac {1}{2} a^3 d x^2+\frac {1}{3} a^3 e x^3+\frac {1}{4} a^3 f x^4+\frac {3}{5} a^2 b c x^5+\frac {1}{2} a^2 b d x^6+\frac {3}{7} a^2 b e x^7+\frac {3}{8} a^2 b f x^8+\frac {1}{3} a b^2 c x^9+\frac {3}{10} a b^2 d x^{10}+\frac {3}{11} a b^2 e x^{11}+\frac {1}{4} a b^2 f x^{12}+\frac {1}{13} b^3 c x^{13}+\frac {1}{14} b^3 d x^{14}+\frac {1}{15} b^3 e x^{15}+\frac {1}{16} b^3 f x^{16} \]

input
Integrate[(c + d*x + e*x^2 + f*x^3)*(a + b*x^4)^3,x]
 
output
a^3*c*x + (a^3*d*x^2)/2 + (a^3*e*x^3)/3 + (a^3*f*x^4)/4 + (3*a^2*b*c*x^5)/ 
5 + (a^2*b*d*x^6)/2 + (3*a^2*b*e*x^7)/7 + (3*a^2*b*f*x^8)/8 + (a*b^2*c*x^9 
)/3 + (3*a*b^2*d*x^10)/10 + (3*a*b^2*e*x^11)/11 + (a*b^2*f*x^12)/4 + (b^3* 
c*x^13)/13 + (b^3*d*x^14)/14 + (b^3*e*x^15)/15 + (b^3*f*x^16)/16
 
3.2.48.3 Rubi [A] (verified)

Time = 0.34 (sec) , antiderivative size = 151, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.120, Rules used = {2017, 2188, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \left (a+b x^4\right )^3 \left (c+d x+e x^2+f x^3\right ) \, dx\)

\(\Big \downarrow \) 2017

\(\displaystyle \int \left (e x^2+d x+c\right ) \left (b x^4+a\right )^3dx+\frac {f \left (a+b x^4\right )^4}{16 b}\)

\(\Big \downarrow \) 2188

\(\displaystyle \int \left (b^3 e x^{14}+b^3 d x^{13}+b^3 c x^{12}+3 a b^2 e x^{10}+3 a b^2 d x^9+3 a b^2 c x^8+3 a^2 b e x^6+3 a^2 b d x^5+3 a^2 b c x^4+a^3 e x^2+a^3 d x+a^3 c\right )dx+\frac {f \left (a+b x^4\right )^4}{16 b}\)

\(\Big \downarrow \) 2009

\(\displaystyle a^3 c x+\frac {1}{2} a^3 d x^2+\frac {1}{3} a^3 e x^3+\frac {3}{5} a^2 b c x^5+\frac {1}{2} a^2 b d x^6+\frac {3}{7} a^2 b e x^7+\frac {1}{3} a b^2 c x^9+\frac {3}{10} a b^2 d x^{10}+\frac {3}{11} a b^2 e x^{11}+\frac {f \left (a+b x^4\right )^4}{16 b}+\frac {1}{13} b^3 c x^{13}+\frac {1}{14} b^3 d x^{14}+\frac {1}{15} b^3 e x^{15}\)

input
Int[(c + d*x + e*x^2 + f*x^3)*(a + b*x^4)^3,x]
 
output
a^3*c*x + (a^3*d*x^2)/2 + (a^3*e*x^3)/3 + (3*a^2*b*c*x^5)/5 + (a^2*b*d*x^6 
)/2 + (3*a^2*b*e*x^7)/7 + (a*b^2*c*x^9)/3 + (3*a*b^2*d*x^10)/10 + (3*a*b^2 
*e*x^11)/11 + (b^3*c*x^13)/13 + (b^3*d*x^14)/14 + (b^3*e*x^15)/15 + (f*(a 
+ b*x^4)^4)/(16*b)
 

3.2.48.3.1 Defintions of rubi rules used

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 2017
Int[(Px_)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[Coeff[Px, x, n - 
 1]*((a + b*x^n)^(p + 1)/(b*n*(p + 1))), x] + Int[(Px - Coeff[Px, x, n - 1] 
*x^(n - 1))*(a + b*x^n)^p, x] /; FreeQ[{a, b}, x] && PolyQ[Px, x] && IGtQ[p 
, 1] && IGtQ[n, 1] && NeQ[Coeff[Px, x, n - 1], 0] && NeQ[Px, Coeff[Px, x, n 
 - 1]*x^(n - 1)] &&  !MatchQ[Px, (Qx_.)*((c_) + (d_.)*x^(m_))^(q_) /; FreeQ 
[{c, d}, x] && PolyQ[Qx, x] && IGtQ[q, 1] && IGtQ[m, 1] && NeQ[Coeff[Qx*(a 
+ b*x^n)^p, x, m - 1], 0] && GtQ[m*q, n*p]]
 

rule 2188
Int[(Pq_)*((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Int[Expand 
Integrand[Pq*(a + b*x + c*x^2)^p, x], x] /; FreeQ[{a, b, c}, x] && PolyQ[Pq 
, x] && IGtQ[p, -2]
 
3.2.48.4 Maple [A] (verified)

Time = 1.50 (sec) , antiderivative size = 151, normalized size of antiderivative = 1.00

method result size
gosper \(a^{3} c x +\frac {1}{2} a^{3} d \,x^{2}+\frac {1}{3} a^{3} e \,x^{3}+\frac {1}{4} f \,a^{3} x^{4}+\frac {3}{5} a^{2} b c \,x^{5}+\frac {1}{2} a^{2} b d \,x^{6}+\frac {3}{7} a^{2} b e \,x^{7}+\frac {3}{8} f \,a^{2} b \,x^{8}+\frac {1}{3} a \,b^{2} c \,x^{9}+\frac {3}{10} a \,b^{2} d \,x^{10}+\frac {3}{11} a \,b^{2} e \,x^{11}+\frac {1}{4} a \,b^{2} f \,x^{12}+\frac {1}{13} b^{3} c \,x^{13}+\frac {1}{14} b^{3} d \,x^{14}+\frac {1}{15} b^{3} e \,x^{15}+\frac {1}{16} b^{3} f \,x^{16}\) \(151\)
default \(a^{3} c x +\frac {1}{2} a^{3} d \,x^{2}+\frac {1}{3} a^{3} e \,x^{3}+\frac {1}{4} f \,a^{3} x^{4}+\frac {3}{5} a^{2} b c \,x^{5}+\frac {1}{2} a^{2} b d \,x^{6}+\frac {3}{7} a^{2} b e \,x^{7}+\frac {3}{8} f \,a^{2} b \,x^{8}+\frac {1}{3} a \,b^{2} c \,x^{9}+\frac {3}{10} a \,b^{2} d \,x^{10}+\frac {3}{11} a \,b^{2} e \,x^{11}+\frac {1}{4} a \,b^{2} f \,x^{12}+\frac {1}{13} b^{3} c \,x^{13}+\frac {1}{14} b^{3} d \,x^{14}+\frac {1}{15} b^{3} e \,x^{15}+\frac {1}{16} b^{3} f \,x^{16}\) \(151\)
norman \(a^{3} c x +\frac {1}{2} a^{3} d \,x^{2}+\frac {1}{3} a^{3} e \,x^{3}+\frac {1}{4} f \,a^{3} x^{4}+\frac {3}{5} a^{2} b c \,x^{5}+\frac {1}{2} a^{2} b d \,x^{6}+\frac {3}{7} a^{2} b e \,x^{7}+\frac {3}{8} f \,a^{2} b \,x^{8}+\frac {1}{3} a \,b^{2} c \,x^{9}+\frac {3}{10} a \,b^{2} d \,x^{10}+\frac {3}{11} a \,b^{2} e \,x^{11}+\frac {1}{4} a \,b^{2} f \,x^{12}+\frac {1}{13} b^{3} c \,x^{13}+\frac {1}{14} b^{3} d \,x^{14}+\frac {1}{15} b^{3} e \,x^{15}+\frac {1}{16} b^{3} f \,x^{16}\) \(151\)
risch \(a^{3} c x +\frac {1}{2} a^{3} d \,x^{2}+\frac {1}{3} a^{3} e \,x^{3}+\frac {1}{4} f \,a^{3} x^{4}+\frac {3}{5} a^{2} b c \,x^{5}+\frac {1}{2} a^{2} b d \,x^{6}+\frac {3}{7} a^{2} b e \,x^{7}+\frac {3}{8} f \,a^{2} b \,x^{8}+\frac {1}{3} a \,b^{2} c \,x^{9}+\frac {3}{10} a \,b^{2} d \,x^{10}+\frac {3}{11} a \,b^{2} e \,x^{11}+\frac {1}{4} a \,b^{2} f \,x^{12}+\frac {1}{13} b^{3} c \,x^{13}+\frac {1}{14} b^{3} d \,x^{14}+\frac {1}{15} b^{3} e \,x^{15}+\frac {1}{16} b^{3} f \,x^{16}\) \(151\)
parallelrisch \(a^{3} c x +\frac {1}{2} a^{3} d \,x^{2}+\frac {1}{3} a^{3} e \,x^{3}+\frac {1}{4} f \,a^{3} x^{4}+\frac {3}{5} a^{2} b c \,x^{5}+\frac {1}{2} a^{2} b d \,x^{6}+\frac {3}{7} a^{2} b e \,x^{7}+\frac {3}{8} f \,a^{2} b \,x^{8}+\frac {1}{3} a \,b^{2} c \,x^{9}+\frac {3}{10} a \,b^{2} d \,x^{10}+\frac {3}{11} a \,b^{2} e \,x^{11}+\frac {1}{4} a \,b^{2} f \,x^{12}+\frac {1}{13} b^{3} c \,x^{13}+\frac {1}{14} b^{3} d \,x^{14}+\frac {1}{15} b^{3} e \,x^{15}+\frac {1}{16} b^{3} f \,x^{16}\) \(151\)

input
int((f*x^3+e*x^2+d*x+c)*(b*x^4+a)^3,x,method=_RETURNVERBOSE)
 
output
a^3*c*x+1/2*a^3*d*x^2+1/3*a^3*e*x^3+1/4*f*a^3*x^4+3/5*a^2*b*c*x^5+1/2*a^2* 
b*d*x^6+3/7*a^2*b*e*x^7+3/8*f*a^2*b*x^8+1/3*a*b^2*c*x^9+3/10*a*b^2*d*x^10+ 
3/11*a*b^2*e*x^11+1/4*a*b^2*f*x^12+1/13*b^3*c*x^13+1/14*b^3*d*x^14+1/15*b^ 
3*e*x^15+1/16*b^3*f*x^16
 
3.2.48.5 Fricas [A] (verification not implemented)

Time = 0.27 (sec) , antiderivative size = 150, normalized size of antiderivative = 0.99 \[ \int \left (c+d x+e x^2+f x^3\right ) \left (a+b x^4\right )^3 \, dx=\frac {1}{16} \, b^{3} f x^{16} + \frac {1}{15} \, b^{3} e x^{15} + \frac {1}{14} \, b^{3} d x^{14} + \frac {1}{13} \, b^{3} c x^{13} + \frac {1}{4} \, a b^{2} f x^{12} + \frac {3}{11} \, a b^{2} e x^{11} + \frac {3}{10} \, a b^{2} d x^{10} + \frac {1}{3} \, a b^{2} c x^{9} + \frac {3}{8} \, a^{2} b f x^{8} + \frac {3}{7} \, a^{2} b e x^{7} + \frac {1}{2} \, a^{2} b d x^{6} + \frac {3}{5} \, a^{2} b c x^{5} + \frac {1}{4} \, a^{3} f x^{4} + \frac {1}{3} \, a^{3} e x^{3} + \frac {1}{2} \, a^{3} d x^{2} + a^{3} c x \]

input
integrate((f*x^3+e*x^2+d*x+c)*(b*x^4+a)^3,x, algorithm="fricas")
 
output
1/16*b^3*f*x^16 + 1/15*b^3*e*x^15 + 1/14*b^3*d*x^14 + 1/13*b^3*c*x^13 + 1/ 
4*a*b^2*f*x^12 + 3/11*a*b^2*e*x^11 + 3/10*a*b^2*d*x^10 + 1/3*a*b^2*c*x^9 + 
 3/8*a^2*b*f*x^8 + 3/7*a^2*b*e*x^7 + 1/2*a^2*b*d*x^6 + 3/5*a^2*b*c*x^5 + 1 
/4*a^3*f*x^4 + 1/3*a^3*e*x^3 + 1/2*a^3*d*x^2 + a^3*c*x
 
3.2.48.6 Sympy [A] (verification not implemented)

Time = 0.03 (sec) , antiderivative size = 180, normalized size of antiderivative = 1.19 \[ \int \left (c+d x+e x^2+f x^3\right ) \left (a+b x^4\right )^3 \, dx=a^{3} c x + \frac {a^{3} d x^{2}}{2} + \frac {a^{3} e x^{3}}{3} + \frac {a^{3} f x^{4}}{4} + \frac {3 a^{2} b c x^{5}}{5} + \frac {a^{2} b d x^{6}}{2} + \frac {3 a^{2} b e x^{7}}{7} + \frac {3 a^{2} b f x^{8}}{8} + \frac {a b^{2} c x^{9}}{3} + \frac {3 a b^{2} d x^{10}}{10} + \frac {3 a b^{2} e x^{11}}{11} + \frac {a b^{2} f x^{12}}{4} + \frac {b^{3} c x^{13}}{13} + \frac {b^{3} d x^{14}}{14} + \frac {b^{3} e x^{15}}{15} + \frac {b^{3} f x^{16}}{16} \]

input
integrate((f*x**3+e*x**2+d*x+c)*(b*x**4+a)**3,x)
 
output
a**3*c*x + a**3*d*x**2/2 + a**3*e*x**3/3 + a**3*f*x**4/4 + 3*a**2*b*c*x**5 
/5 + a**2*b*d*x**6/2 + 3*a**2*b*e*x**7/7 + 3*a**2*b*f*x**8/8 + a*b**2*c*x* 
*9/3 + 3*a*b**2*d*x**10/10 + 3*a*b**2*e*x**11/11 + a*b**2*f*x**12/4 + b**3 
*c*x**13/13 + b**3*d*x**14/14 + b**3*e*x**15/15 + b**3*f*x**16/16
 
3.2.48.7 Maxima [A] (verification not implemented)

Time = 0.20 (sec) , antiderivative size = 150, normalized size of antiderivative = 0.99 \[ \int \left (c+d x+e x^2+f x^3\right ) \left (a+b x^4\right )^3 \, dx=\frac {1}{16} \, b^{3} f x^{16} + \frac {1}{15} \, b^{3} e x^{15} + \frac {1}{14} \, b^{3} d x^{14} + \frac {1}{13} \, b^{3} c x^{13} + \frac {1}{4} \, a b^{2} f x^{12} + \frac {3}{11} \, a b^{2} e x^{11} + \frac {3}{10} \, a b^{2} d x^{10} + \frac {1}{3} \, a b^{2} c x^{9} + \frac {3}{8} \, a^{2} b f x^{8} + \frac {3}{7} \, a^{2} b e x^{7} + \frac {1}{2} \, a^{2} b d x^{6} + \frac {3}{5} \, a^{2} b c x^{5} + \frac {1}{4} \, a^{3} f x^{4} + \frac {1}{3} \, a^{3} e x^{3} + \frac {1}{2} \, a^{3} d x^{2} + a^{3} c x \]

input
integrate((f*x^3+e*x^2+d*x+c)*(b*x^4+a)^3,x, algorithm="maxima")
 
output
1/16*b^3*f*x^16 + 1/15*b^3*e*x^15 + 1/14*b^3*d*x^14 + 1/13*b^3*c*x^13 + 1/ 
4*a*b^2*f*x^12 + 3/11*a*b^2*e*x^11 + 3/10*a*b^2*d*x^10 + 1/3*a*b^2*c*x^9 + 
 3/8*a^2*b*f*x^8 + 3/7*a^2*b*e*x^7 + 1/2*a^2*b*d*x^6 + 3/5*a^2*b*c*x^5 + 1 
/4*a^3*f*x^4 + 1/3*a^3*e*x^3 + 1/2*a^3*d*x^2 + a^3*c*x
 
3.2.48.8 Giac [A] (verification not implemented)

Time = 0.28 (sec) , antiderivative size = 150, normalized size of antiderivative = 0.99 \[ \int \left (c+d x+e x^2+f x^3\right ) \left (a+b x^4\right )^3 \, dx=\frac {1}{16} \, b^{3} f x^{16} + \frac {1}{15} \, b^{3} e x^{15} + \frac {1}{14} \, b^{3} d x^{14} + \frac {1}{13} \, b^{3} c x^{13} + \frac {1}{4} \, a b^{2} f x^{12} + \frac {3}{11} \, a b^{2} e x^{11} + \frac {3}{10} \, a b^{2} d x^{10} + \frac {1}{3} \, a b^{2} c x^{9} + \frac {3}{8} \, a^{2} b f x^{8} + \frac {3}{7} \, a^{2} b e x^{7} + \frac {1}{2} \, a^{2} b d x^{6} + \frac {3}{5} \, a^{2} b c x^{5} + \frac {1}{4} \, a^{3} f x^{4} + \frac {1}{3} \, a^{3} e x^{3} + \frac {1}{2} \, a^{3} d x^{2} + a^{3} c x \]

input
integrate((f*x^3+e*x^2+d*x+c)*(b*x^4+a)^3,x, algorithm="giac")
 
output
1/16*b^3*f*x^16 + 1/15*b^3*e*x^15 + 1/14*b^3*d*x^14 + 1/13*b^3*c*x^13 + 1/ 
4*a*b^2*f*x^12 + 3/11*a*b^2*e*x^11 + 3/10*a*b^2*d*x^10 + 1/3*a*b^2*c*x^9 + 
 3/8*a^2*b*f*x^8 + 3/7*a^2*b*e*x^7 + 1/2*a^2*b*d*x^6 + 3/5*a^2*b*c*x^5 + 1 
/4*a^3*f*x^4 + 1/3*a^3*e*x^3 + 1/2*a^3*d*x^2 + a^3*c*x
 
3.2.48.9 Mupad [B] (verification not implemented)

Time = 9.26 (sec) , antiderivative size = 150, normalized size of antiderivative = 0.99 \[ \int \left (c+d x+e x^2+f x^3\right ) \left (a+b x^4\right )^3 \, dx=\frac {f\,a^3\,x^4}{4}+\frac {e\,a^3\,x^3}{3}+\frac {d\,a^3\,x^2}{2}+c\,a^3\,x+\frac {3\,f\,a^2\,b\,x^8}{8}+\frac {3\,e\,a^2\,b\,x^7}{7}+\frac {d\,a^2\,b\,x^6}{2}+\frac {3\,c\,a^2\,b\,x^5}{5}+\frac {f\,a\,b^2\,x^{12}}{4}+\frac {3\,e\,a\,b^2\,x^{11}}{11}+\frac {3\,d\,a\,b^2\,x^{10}}{10}+\frac {c\,a\,b^2\,x^9}{3}+\frac {f\,b^3\,x^{16}}{16}+\frac {e\,b^3\,x^{15}}{15}+\frac {d\,b^3\,x^{14}}{14}+\frac {c\,b^3\,x^{13}}{13} \]

input
int((a + b*x^4)^3*(c + d*x + e*x^2 + f*x^3),x)
 
output
(a^3*d*x^2)/2 + (b^3*c*x^13)/13 + (a^3*e*x^3)/3 + (b^3*d*x^14)/14 + (a^3*f 
*x^4)/4 + (b^3*e*x^15)/15 + (b^3*f*x^16)/16 + a^3*c*x + (3*a^2*b*c*x^5)/5 
+ (a*b^2*c*x^9)/3 + (a^2*b*d*x^6)/2 + (3*a*b^2*d*x^10)/10 + (3*a^2*b*e*x^7 
)/7 + (3*a*b^2*e*x^11)/11 + (3*a^2*b*f*x^8)/8 + (a*b^2*f*x^12)/4